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Abstract 

A simple d-dimensional lattice model is proposed, incorporating some degree 
of frustration and thus capable of describing some aspects of molecular orien- 
tation in covalently bound molecular solids. For d = 2 the model is shown to 
be equivalent to the standard two-dimensional Ising model, while for d > 2 it 
describes a peculiar transition from an isotropic high temperature phase to a 
low-dimensional anisotropic low temperature state. A general mean field analysis 
is presented and compared to some exact limit properties. 
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Most molecular liquids retain their molecular structure even in the solid phase, where 
some long range order usually shows up as a consequence of inter-molecular interaction. 
However in the solid the orientational order of the molecules may change according to the 
thermodynamic conditions giving rise to quite rich phase diagrams as recently observed for 
hydrogen under high pressure [1]. Orientational ordering is responsible for several phase 
transitions occurring even in the liquid phase (liquid crystals) and is thought to be related 
to the metal-insulator transition of liquid iodine [2] . 

Molecular ordering may be classified into two different classes according to the nature 
of the interaction. The usual interaction has a quadrupole nature, is weak and gives 
rise to the observed three-dimensional ordering of most molecular Van der Waals solids. 
Conversely in some molecular solids the interaction has a covalent main component and is 
characterized by some level of frustration since the coordination number for the covalent 
bond is quite low. Each molecule must choose a few partners and cannot accept any 
further proposal. The lower is the allowed coordination number, the higher the frustration 
which gives rise to a low-dimensional structure as observed in polymers (one-dimensional) 
or in solid iodine [3,4] (two-dimensional). Actually some one-dimensional zigzag chains 
have also been reported [5] in iodine where the covalent nature of the interaction is out 
of doubt [6,7]. Moreover we expect that a covalent interaction should show up for all 
the molecular solids under high pressure, as the inter-molecular distance approaches the 
intra-molecular length, provided that some important structural transition does not occur 
first (as dissociation). 

While the thermodynamic behaviour of the quadrupole driven ordering can be an- 
alyzed in terms of 0(3) symmetric vectorial models, in this letter we introduce a sim- 
ple frustrated lattice model which is capable of describing some aspects of a covalently 
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bound molecular solid. More precisely the model describes the transition from an high- 
temperature (or weakly interacting) fully isotropic disordered system, to a low-temperature 
(or strongly interacting) anisotropic low dimensional broken-symmetry phase. As a conse- 
quence of frustration the breaking of symmetry is accompanied by a sort of decomposition 
of the system in low-dimensional almost independent parts, as observed in solid iodine. 
Such remarkable behaviour requires a space dimension d > 2, while for d = 2 the model 
is shown to be equivalent to the exactly solvable two-dimensional Ising model [8]. Then, 
apart from its original motivation, the model seems to be of interest by itself as an uncon- 
ventional generalization of the standard Ising model to higher dimensions. 

Let us consider a <i-dimensional hypercubic lattice, with a randomly oriented linear 
molecule at each site. The molecules are supposed to be symmetric with respect to their 
centre of mass which is fixed at the lattice site. The description is simplified by allowing 
only a discrete number of space orientations for each molecule: we assume that each of 
them must point towards one of its 2d first neighbour sites. This choice can be justified by 
the existence of covalent interactions along preferred axes. In this quantized version each 
molecule has d different states corresponding to molecular orientation along the hypercube 
axes (molecules are supposed to be symmetric). Finally, each couple of first neighbour 
molecules, when pointing the one towards the other, are assumed to gain a bonding energy 
for their directional covalent bond (they touch each other). As shown in figure 1 for d = 2, 
bonding in a direction excludes any possible bond along the other (d — 1) directions. The 
coordination number is 2 for any value of d, and the frustration increases with increasing 
d. 

According to such description we introduce a variable versor w r for each of the N 
sites r of the lattice, with w r G {xi,X2, ' • ' x d\ pointing towards one of the d hypercube 
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The partition function follows 

Z = e S = ^ ^ r ' a **'* a 5 * r+ * a &a (!) 

{w} {w} 

where {w} indicates a sum over all configurations, a runs from 1 to d, and the lattice 
spacing is set to unity. The inverse temperature ft, in units of binding energy, can be 
negative for a repulsive model, but is assumed positive in the present molecular context. 

We first show the equivalence for d = 2 between molecular and Ising model. Let us 
add for each couple of sites a generic field h(a) according to the space direction a of the 
lattice link which joins the sites. The partition function reads 

4/3^ (^ r i« • S w r+&a x a + h(a)w r + h(a)w r+£a ) 

z h = J2 eSh = J2 e r,Q ( 2 ) 

{w} {w} 



Introducing a matrix notation 



^ = X> 5h = X> Er - L(r,a) ( 3 ) 

{w} {w} 

the Lagrangian density L follows as 

L(r, a) = wlM a (P, h)w r+£a (4) 

Here the canonical rf-dimensional column vector representation of R d is employed with 
xi = (1,0,0...), i-2 = (0,1,0...), etc. The dx d matrix M a does not depend on the 
configurations of the system, and entirely characterizes it. Now provided that ^ a h(a) = 
0, Sh does not depend on h and Sh = S. For d = 2 such condition is satisfied by the field 
h(l) = hfa - x 2 ), h(2) = -h(l). The matrix M a follows 

Mi /4 W + 2ft) o \, M2 = (-m \ 



-8/3/i / ' 1 V 4/3(1 + 2/i) 



Then for h = -1/4, M x = M 2 , and L reads 



L(r, a) = /3 + u)J 



t f -0 



) 



a 



(6) 



Identifying the two-dimensional versors w with spin variables, apart from an inessential 
factor, Z reduces to the partition function of a two-dimensional Ising model 



and is exactly solvable. For (3 — > +oo a ground state is approached with all the molecules 
aligned along the same direction, and formation of one-dimensional polymeric chains (figure 
2a); for (3 — > — oo the repulsive model approaches a zero-energy (no bonds) ground state 
analogous to the antiferromagnetic configuration of the Ising model (figure 2b). 

For d > 3 the analogy with the Ising model breaks down, and this is evident from 
a simple analysis of the ground state configuration. Due to frustration the model has an 
infinitely degenerate ground state in the thermodynamic limit N — > oo. For instance in the 
case d = 3, the minimum energy is obtained by aligning all the molecules along a common 
direction, as for d = 2. However the ground state configuration is not unique: the number 
of molecular bonds does not change if we rotate together all the molecules belonging to 
an entire layer which is parallel to the original direction of alignment. As a consequence 
of frustration the total degeneration is 3(2^ /3 )), and the system is expected to behave 
like a glass for the large energy threshold which separate each minimum from the other. 
The phase diagram is expected to be quite rich, with at least a transition point between 
the high temperature disordered phase and an ordered broken-symmetry low temperature 
phase. 



Z = e? Nd ■ Z 



Ising 



(7) 
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Some analytical result can be obtained in Mean-Field (MF) approximation: neglecting 
second order fluctuation terms 



$w r x a ' ^w r+S: x a ~ $w r x a ' A a + ^w r +xo,x a ' A a (8) 



where A Q = (8n, r $ a ), and ^ Q A a = 1 (with the constraint < A a < 1). The partition 
function factorizes as 

.i 



Zmf = e 47V/3 S„ A c 
and the free energy follows 



e 

a 



JV 

8/3A 



0) 



F«F = log Z MF = 4 £ - i log ^- \ . 

The derivative with respect to A M yields, for the stationary points 



(10) 



e 



8/3A M 



which satisfies the condition ^ Q A Q = 1. 

In the high temperature limit (3 — > eq. (11) admits the unique solution A^ = 1/d 
which reflects the complete random orientation of molecules. In the opposite limit (3 — > oo, 
apart from such solution, eq. (11) is satisfied by the broken-symmetry field A^ = 1, 
A a = for a ^ fx, which obviously corresponds to a minimum for Fmf- Then at a critical 
point (3 = (3 C the high temperature solution must become unstable towards a multivalued 
minimum configuration. The Hessian matrix is easily evaluated at the stationary points 
employing eq. (11) 

H ^ = l^x^k = ^ (1 ~ 8/3A ^ } + 8/3A ^ (12) 



For (3 < P c , inserting A^ = 1/d, the eigenvalue problem det\H^ — \5^\ = yields 



Thus the Hessian matrix is positive defined if and only if A = (1 — 8/3 /d) > 0. Beyond the 
critical point (3 = f3 c = (d/8) the solution A^ = 1/d is not a minimum, and a multivalued 
minimum configuration shows up. Such result obviously agrees with the MF prediction for 
the Ising model I3i sing = l/(2rf), only for the special dimension d = 2. For d > 2 we observe 
an increase of (3 C with d, to be compared to the opposite trend shown by the Ising model. 
Such behaviour may be interpreted in terms of the low dimensionality of the ordered phase. 
Due to frustration the ordering may only occur on a low dimensional scale: for instance 
in three dimensions each layer has an independent internal ordering. Thus we expect a 
larger j3 c for d > 2 since the increasing of d only introduces larger fluctuations, with each 
molecule having (d — 2) allowed out-of-plane orientations. For d = 3 the low temperature 
phase can be regarded as a quenched disordered superposition of layers which are internally 
ordered along different in-plane directions. Thus as a consequence of frustration the system 
shows a two-dimensional character below the critical point while behaving as truly three- 
dimensional in the high temperature domain. 

We should wonder if the MF approach is reliable and how much this description is 
affected by the MF approximation. In fact MF associates the existence of long range order 
with some broken-symmetry anisotropic field. On the other hand frustration enforces the 
long range order to be anisotropic in space: this is exact since we have shown that, at 
variance with usual magnetic models, even the ground state does not have an isotropic 
long range order. Long range order may only occur inside low-dimensional sub-sets of the 
lattice. Such sub-sets may even be one-dimensional, as infinite polymeric chains of aligned 




(13) 
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molecules. For d = 2 the existence of such chains is equivalent to the existence of two- 
dimensional long range order; for d > 2 the existence of infinitely long one-dimensional 
chains does not imply the existence of a two-dimensional order, since the chains could be 
randomly oriented along all the lattice axes. In such sense long range order does not neces- 
sarily imply the breaking of isotropy. Any configuration characterized by one-dimensional 
long range order would have an higher energy compared to the ground state, and could be 
stable only in an intermediate range of temperature. We expect that decreasing the tem- 
perature the generic d-dimensional system could reach first a one-dimensional order and 
then the full two-dimensional order (i.e. order inside each two-dimensional layer). This 
already corresponds to a minimum of the configurational energy, since nothing is gained 
by extending the ordering to higher dimensional sub-sets. The possible existence of an 
intermediate isotropic phase characterized by one-dimensional long range order is an open 
question which cannot be addressed by the present MF approach. The critical tempera- 
ture should be higher as compared to full ordering since a free one-dimensional chain of 
molecules, each having d allowed orientational states, always shows long range order. By 
an argument similar to that employed below eq. (4), such one-dimensional chain can be 
mapped on a one-dimensional magnetic model in presence of a non-zero magnetic field, 
and is characterized by a non-zero net magnetization for any value of (3. 

Such speculation is corroborated by several estimates of (3 C obtained for d = 3 by a 
modified Migdal-Kadanoff renormalization group approach [9] or by employing an approx- 
imate duality relation [10]. In both cases no anisotropic field is required, and the transition 
points are detected through the singular behaviour of the free energy at the onset of long 
range order. Details about such techniques will be published elsewhere. We may anticipate 
the prediction (3 C xs 0.15 for d = 3, to be compared to the MF result (3mf = 3/8 = 0.375. 
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Since any improving of the MF approach would eventually yield an increase of (3m f as a 
consequence of fluctuations [10], we are tempted to attribute the different estimates for j3 c 
to two different transitions. Of course further work is required in order to establish the 
possible existence of such intermediate phase. 

In summary a new non-trivial extension of the Ising model has been proposed for 
d > 2, incorporating frustration effects and thus reproducing some aspects of ordering in 
covalent molecular systems. The model undergoes an order-disorder transition from an 
high temperature isotropic state to a low temperature anisotropic low-dimensional phase. 
A MF analysis is presented for a generic d > 2, and the results have been discussed 
in comparison to some exact limit properties of the <i-dimensional model. The possible 
existence of an intermediate isotropic phase, characterized by one-dimensional long range 
order, is proposed as an open question for further work on the model. 
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Figure Captions 

Fig.l An allowed configuration for d = 2. 

Fig. 2 Ground state configurations for the two-dimensional attractive (a) and repulsive 
models. 
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